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ABSTRACT 

There is growing evidence from both spectral and timing properties that there is a 
truncated inner accretion disc in low mass accretion rate Galactic black hole systems. 
The detection of extremely smeared relativistic iron lines in some of these systems is 
the only current piece of evidence which conflicts with this geometrical interpretation 
of the low/hard state. Here we show that the line width in the BeppoSAX data of a 
bright low/hard state of the transient black hole XTE J 1650-500 is indeed consistent 
with extreme relativistic effects. However, the relativistic smearing can be significantly 
reduced if there is also resonance iron K line absorption from an outflowing disc wind. 
The iron line smearing is then completely compatible with a truncated disc, so gives 
no information on the black hole spin. 

Key words: accretion, accretion discs - black hole physics - X-rays: binaries - 
X-rays: individual: XTE J1650-500 



1 INTRODUCTION 

The existence of a minimum stable orbit for material around 
a black hole is a key prediction of general relativity in the 
strong field limit. For Schwarzschild black holes this is at 
Rms = 6i? 9 (where R g = GM/c 2 ), while for maximally ro- 
tating Kerr black holes (spin parameter a — 0.998) this re- 
duces to 1.23-Rg. Bright accretion flows give a way to obser- 
vationally test such ideas as well as to directly estimate the 
spin of the black hole. Optically thick material in Keplerian 
orbits should emit a quasi-thermal spectrum with maximum 
temperature emission from the largest luminosity /smallest 
area regions i.e. set by the last stable orbit (Shakura & Sun- 
yaev 1973). As these models are thermal, there is a clear 
prediction that their maximum temperature, T max , should 
increase with total disc luminosity as Ldisc oc T^, ax , if the 
emission is from a constant area i.e. if there is a constant in- 
ner disc radius. Strong observational support for this comes 
from the Galactic black hole (GBH) binary systems. These 
can show disc dominated spectra (often termed the high/soft 
state) where Ldisc oc T^ ax over large changes in luminosity 
(Ebisawa et al. 1993; Kubota, Makishima & Ebisawa 2001; 
Kubota & Makishima 2004) . The fixed radius inferred from 
the proportionality constant of this relation is consistent 
with a Keplerian disc (stress-free inner boundary condition) 
in a low-to-moderate spin spacetime (a < 0.8) for all GBH 
for which this experiment can be performed (Gierlihski & 
Done 2004; Li et al. 2005; Shafee et al. 2005; Davis, Done 



& Blaes 2005). This ties in with the spins expected from 
models of stellar collapse. These predict a maximum birth 
spin of a ~ 0.8 from the combination of slow rotation of the 
pre-supernovae core and angular momentum losses through 
gravitational wave radiation during collapse (e.g. Gammie, 
Shapiro & McKinney 2004), and there is not enough mass 
transfer in the binary evolution to cause appreciable spin up 
from accretion (King & Kolb 1996). 

However, there is growing evidence that the accretion 
disc inner radius can also be larger than 7? ms , due to the 
inner disc evaporating into hot, geometrically thick accre- 
tion flow. These models were developed in response to ob- 
servations of GBH at low luminosity (generally termed the 
low/hard state), where the expected accretion disc emission 
is much weaker and at lower temperature than predicted 
by an optically thick disc extending down to the last stable 
orbit (e.g Esin et al. 2001; Frontera et al. 2001). Instead, 
the X-ray spectra peak at ~100 keV (see e.g. the review 
by Zdziarski & Gierlihski 2004), requiring that the contin- 
uum is formed in hot, optically thin plasma, which is not 
strongly cooled by Compton scattering of seed photons (e.g. 
Zdziarski & Gierlihski 2004). These cooling limits are very 
difficult (Barrio, Done & Nayakshin 2003; though not im- 
possible e.g. Beloborodov 1999) to reconcile with a cool disc 
extending down to the minimum stable orbit. Other obser- 
vational advantages of truncated disc configuration is that 
it can provide a geometric explanation for the major X-ray 
spectral state transitions seen in GBH by a physical change 
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in the accretion flow structure. At low luminosities the in- 
ner cool disc is replaced by a hot accretion flow. The tran- 
sition radius between the two phases is far from the hole, 
few photons from the disc illuminate the hot flow and the 
spectrum is hard. As the transition radius moves inwards 
the disc penetrates further into the hot flow, so more seed 
photons cool the flow and its spectrum softens until the disc 
extends down to the last stable orbit (Esin et al. 1997). This 
changing transition radius also gives a changing size scale to 
associate with the changing characteristic frequencies seen 
in the variability power spectra (see e.g. the review by Done 
2001). 

Theoretically, such hot inner flows are generically pro- 
duced by a range of analytic approximations to the accretion 
flow equations. In these, the protons are heated by grav- 
ity, while the electrons gain energy only through Coulomb 
collisions and cool predominantly via Compton scattering 
(Shapiro, Lightman & Eardley 1975), advection (Narayan 
& Yi 1995), outflows (Blandford & Begelman 1999), con- 
vection (Abramowicz & Igumenshchev 2001) or jets (Fal- 
cke, Kording & Markoff 2004). Numerical simulations of 
the accretion flow at low densities, including the self con- 
sistently generated magnetic dynamo viscosity (Balbus & 
Hawley 1991), appear to show some form of inner hot flow, 
though its properties are more complex than any of the ap- 
proximate analytic solutions above (Hawley & Balbus 2002) . 

Another argument for the truncated disc comes from 
the quasi-periodic oscillations (QPOs) observed in the power 
density spectra of black hole binaries. During the transition 
from the low/hard to the high/soft spectral state their fre- 
quency typically increases (e.g. Cui et al. 1999; Rossi et al. 
2004) indicating a decreasing size-scale in the accretion flow. 
If the QPO frequency is directly related to the transition 
radius between the cold disc and the hot inner flow (Stella 
& Vietri 1998; Titarchuk, Osherovich & Kuznetsov 1999; 
Psaltis & Norman 2000) then its behaviour is consistent with 
a decreasing truncation radius as the source approaches the 
soft state. 

Clearly, the point at which the cool disc truncates is of 
great importance to our understanding of the nature of the 
accretion flow, and hence whether we can use observations 
of compact objects to test strong gravity. Fortunately there 
is an independent constraint on this from the shape of the 
fluorescent iron line produced by X-ray illumination of the 
cool disc. Special and general relativistic effects are both 
important as the disc material is in very high velocity (Kep- 
lerian) orbits in a strong gravitational field. All these effects 
get stronger as the disc extends closer to the black hole, 
so the width of the line gives a measure of the inner edge 
of the cool fluorescing material (Fabian et al. 1989; 2000). 
Thus for the accretion flow model above in which there is a 
variable truncation radius between a disc and inner hot flow, 
the line width should be correlated with spectral hardness. 
Hard spectra (generally seen at low L/L E(M ) should have 
narrow lines, softer spectra should have broader lines. There 
are indications that this correlation is seen both in AGN 
(Lubihski & Zdziarski 2001) and in GBH (Zycki, Done & 
Smith 1998; Churazov, Gilfanov & Revnivtsev 2000; Ibrag- 
imov et al. 2005). 

However, some of the more recent results from the line 
shape seen in the GBH do not fit easily into this picture. 
There are reported lines in the low/hard state (i.e hard X-ray 



spectra, with QPO frequencies which are not at their maxi- 
mum values) which are so broad as to require a disc around 
a maximally spinning black hole, illuminated by highly cen- 
trally concentrated radiation field, perhaps implying direct 
tapping of the spin energy of the black hole or at least sub- 
stantial continuous stress across the last stable orbit (Mini- 
utti, Fabian & Miller 2004; Miller et al. 2002; Miller et al. 
2004a). Plainly this would rule out the inner hot flow model 
if there is no alternative explanation for the line width. The 
clearest example of this discrepancy is the BeppoSAX data 
of the bright low/hard state of XTE J1650-500. The line 
width in this spectrum is also in conflict with the moder- 
ate spin derived from the luminosity-temperature relation 
obtained when this object shows disc dominated spectra 
(Gierlinski & Done 2004). Thus the extent of the disc in- 
ferred from the iron line width in XTE J1650-500 seems at 
odds with the broad band spectral shape, the QPOs, and 
the disc continuum emission. 

Here we re-examine the BeppoSAX data on XTE 
J1650-500 and show that while these data are indeed consis- 
tent with an extreme relativistic line, they can be fit equally 
well (or better) using a model of a 'normal' relativistic line 
(i.e. from a truncated disc with standard gravitational en- 
ergy release), together with blueshifted highly ionized nar- 
row absorption features (predominantly resonance lines of 
He- and H-like iron). These absorption features distort the 
blue edge of the line, especially where the reflection spec- 
trum is ionized and/or from high inclination material. If the 
presence of such material is confirmed then this removes the 
only observational evidence against a truncated disc inter- 
pretation of the low/hard state in GBH. 



2 THE DATA 

We extract the XTE J1650-500 BeppoSAX archival data 
taken on 2001 September 11 for both the Medium-Energy 
Concentrator Spectrometer (MECS) and Phoswich Detector 
System (PDS) instruments, using the standard (September 
1997) background, auxiliary response (MECS) and redis- 
tribution files (MECS and PDS). We use the same energy 
ranges as Miniutti et al. (2004), i.e. 2.5-10 keV (MECS) and 
13-200 keV (PDS), with 2 per cent systematic errors added 
in quadrature to the PDS data. We checked the publicly 
available Crab PDS data and obtain reduced xV" = 0.9 
for a single power law of energy index 1.13, showing that 
this prescription is an adequate description of the response 
uncertainties. 

We also extract the RXTE Proportional Counter Ar- 
ray (PCA) and High-Energy X-ray Timing Experiment 
(HEXTE) spectra from the entire outburst. We calculate 
the X-ray colours and estimate the bolometric luminosity 
following the method of Done & Gierlinski (2003). We used 
the distance of 4 kpc (Tomsick et al. 2003) and mass of 7.3 
M (Orosz et al. 2004) for the luminosity estimate. We also 
extracted the power-density spectra from PCA observations 
in the 1/128 to 64 Hz frequency band and fitted them by a 
multiple Lorentzian model. 

We use XSPEC version 11.3.1 (Arnaud 1996) for spectral 
fitting. The error of each model parameter is given for a 90 
per cent confidence interval. We fix the absorption column 



© 2002 RAS, MNRAS 000, 1-11 



Truncated disc versus extremely broad iron line 3 



o. i 



0.01 
10 





- 


* 


.■■ 


1 1 1 1 1 

% 
1 

/ 

.1 


- * I* 




1 1 1 1 : 

■'—4 i 

1 , 1 " , " " 1 , 1 : 



1 =- 



2 

7) 1.5 

I 1 

o 

u 0.5 



52160 



52180 



52200 
Time (MJD) 



52220 



52240 



Figure 1. The top panel shows the evolution of the total, absorption corrected, bolometric luminosity (as a function of Eddington) for 
all the RXTE pointed observations of the outburst of XTE J1650-500. The total flux is derived from fitting the PCA and HEXTE data 
with a model of an accretion disc plus its Comptonized emission (together with a broad line and smeared edge to mimic reflection). The 
middle panel shows the low frequency QPO derived from fitting the broad-band power spectra with multiple Lorentzian components, 
while the bottom panel shows the soft (triangles, red in colour) and hard (black squares) X-ray colours. Both QPO and colours are still 
evolving smoothly at the point where the BeppoSAX data were taken (vertical line). In the context of a truncated disc model for the 
low/hard state, this implies that the disc has not yet reached the last stable orbit. 



at N H = 5 x 10 21 cm" 2 in all fits (Miller et al. 2002; Miniutti 
et al. 2004). 



3 TIMING AND SPECTRAL EVOLUTION 
DURING THE OUTBURST 

Fig. 1 shows the estimated bolometric luminosity, the fre- 
quency of the strong low-frequency QPO and X-ray colours 
obtained from the RXTE. The vertical line marks the time 
of the BeppoSAX observation analysed in this paper. The 
QPO frequency during the BeppoSAX observation is plainly 
not the maximum seen from this source. Any association of 
the QPO with the inner radius of the cool disc clearly re- 
quires that the disc does not extend down to the last stable 
orbit, but truncates above this point. The X-ray colours show 
that the BeppoSAX data were taken while the source is still 
smoothly softening, in a position which is associated with 
the softest low/hard states. Given that the spectrum soft- 
ens further after the BeppoSAX observation, models of the 
broad band spectra with a moving inner disc predict that 
the disc truncates above the last stable orbit. Thus both 
the QPOs and broad band spectral shape indicate that the 
innermost parts of a cool disc are not present. This is in 
conflict with the extreme iron line smearing in these same 
data which requires that the disc extends down to the last 
stable orbit of a high spin black hole. 

© 2002 RAS, MNRAS 000, 1-11 



4 MECS 2-10 KEV DATA: 45° INCLINATION 
4.1 Simple models 

With the MECS data we can reproduce the extreme rela- 
tivistic broadening of Miniutti et al. (2004). We use their 
phenomenological model of a smeared edge and iron line 
with extreme Kerr spacetime distortions (the LAOR model 
in xspec) on a continuum consisting of a disc spectrum and 
power law. This gives an adequate fit (v 2 = 70.9/65) to the 
data, and the parameters of the strong iron line (equivalent 
width of I90I20 e V) do indeed point to a very small inner 
disc radius of n n = Ri n /R g — 2.20lg'ig with extreme central 
concentration of the illumination, oc r~ q where q = 4.3^q'3. 
A similarly extreme line is seen using the same model in the 
XMM-Newton data taken a few days after the BeppoSAX 
observation (Miller et al. 2002). 

The XSPEC LAOR code is written for the specific case 
of calculating the relativistic effects on a line which can 
be approximated as an initial delta function. This is iden- 
tical to the transfer function, so relativistic effects on an 
arbitrary initial spectrum can be calculated by simply con- 
volving the initial spectrum with this transfer function. We 
recoded the LAOR model into a convolution model (here- 
after called CONLINe) in xspec and use this to smear a full 
reflected spectrum (the pexriv model in XSPEC with addi- 
tional narrow gaussian line) to replace the phenomenological 
smeared edge/line features in the previous fit. Convolution 
means that the spectrum outside of the observed energy 
band can be important so we always extend the energy range 
used by xspec to 0.2-50 keV. 

We follow Miller et al. (2002) and Miniutti et al. (2004) 
and fix the inclination angle to 45° and assume solar abun- 
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dances. We constrain the radial emissivity index to be be- 
tween 3 and 5 as suggested by these previous fits. This gives 
X 2 = 75.7/65, i.e. slightly worse than the phenomenological 
approximation for reflected continuum emission. The best- 
fitting model is reflection dominated, with Q/2tt — oo, i.e. 
for a model in which the intrinsic continuum is hidden from 
view. However, this is poorly constrained, with the error 
range allowing Q./2n > 0.7. By contrast the line param- 
eters are fairly well constrained, with equivalent width of 
120 ± 25 eV, and inner disc radius of n n = 2.2+g'.i2 an d 
q — 3. 7^0 3 ■ Thus the line parameters do change with the 
addition of a reflected continuum, most significantly in re- 
ducing the equivalent width, although the smearing is also 
affected, with the illumination being not quite so centrally 
concentrated. Detailed results on the line equivalent width 
(and to a lesser extent the shape) derived without includ- 
ing a reflected continuum should be treated with caution. 
Nonetheless, the small disc inner radius is very similar to 
that in the phenomenological fits; there is a very broad fea- 
ture in the data between 4-7 keV as shown in Fig. 2a. 

4.2 Complex reflection 

The reflected component from the pexriv fits in the previ- 
ous section is strongly ionized. This means that the intrinsic 
iron line emission from the disc is not narrow. It should in- 
stead be strongly Comptonized, and this Comptonization 
should also affect the edge in the reflected continuum (Ross, 
Fabian & Young 1999). The neglect of Compton scattering 
below 12 keV in pexriv means that this code has an edge 
which which is artificially deep and sharp for a highly ion- 
ized disc. We replace the pexriv reflected continuum and 
narrow gaussian line with reflected spectra of D. Ballantyne 
from a constant density ionized disc (hereafter CDID; Bal- 
lantyne, Iwasawa & Fabian 2001, based on Ross & Fabian 
1993) which are publicly available as xspec table models. 
These calculate the reflected spectrum from an X-ray illumi- 
nated slab of constant density, including Compton up- and 
down-scattering below 12 keV, with self consistent iron (and 
other elements) line calculation for solar abundance mate- 
rial. These give \ 2 = 93.7/67 for a highly smeared, ionized 
slab, somewhat worse than the (unphysical) fits above, but 
with two fewer free parameters as the line energy and width 
are now calculated self-consistently with the reflected con- 
tinuum. The reflection emission is highly ionized, lg(£/erg 
cm s _1 ) — 3.79to'o8) an d again the spectrum is consistent 
with pure reflected continuum, although the errors are wide 
and only restrict £1/2% > 0.6. 

The best fit again picks out extreme relativistic smear- 
ing, with q = 51"° 2 an d fin = 3.0^g'4- However, the signif- 
icance of detection of relativistic smearing is reduced with 
respect to the reflection models which do not include Comp- 
tonization of the line and edge features. The pexriv based 
reflection models give A\ 2 ~ 190 for removing all smear- 
ing from the model, while the CDID reflection models give 
A X 2 ~ 13. 

Figs. 2(b) and (c) show these fits with extreme relativis- 
tic smearing, and no relativistic smearing, respectively. It is 
clear that the intrinsic reflected emission alone - mainly the 
edge in the reflected continuum rather than the line itself, 
as there is very little line emission at these high ionizations 
- can match most of the broad feature in the data. This can 



have a strong impact on the detection significance of rela- 
tivistic smearing and some impact on the detailed parame- 
ters. Plainly it is important to use an accurate description of 
ionized reflection. However, there is considerable uncertainty 
in these models as they depend on the (unknown) vertical 
and radial ionization structure of the disc. The CDID mod- 
els assume that the reflecting material can be approximated 
by a slab of constant (vertical and radial) density. Yet a 
disc is in hydrostatic equilibrium so its structure responds 
to the incident radiation (Nayakshin, Kazanas & Kallman 
2000; Ballantyne, Ross & Fabian 2001). The upper layers 
are heated to the Compton temperature of the illuminating 
flux, so they expand, becoming less dense and more ionized. 
For hard power law illumination the Compton temperature 
is large, so the upper layers are completely ionized, forming 
a skin on top of the disc, but for softer illuminating spec- 
tra (such as that of the data used here) the upper layers of 
the disc are still visible so the difference between the hy- 
drostatic and constant ionization reflection models is much 
smaller (Done & Nayakshin 2001). As well as this vertical 
structure, there should also be a radial dependance of the 
ionization structure of the disc, and the full reflected emis- 
sion should be the emissivity weighted integral of all of these 
different radii regions with vertical ionization structure. 

We use the xion code (Nayakshin et al. 2000) to model 
reflection from a disc of material in hydrostatic equilibrium 
from power law illumination, i.e. including both radial and 
vertical ionization gradients. This can be used to produce the 
total emission (direct plus reflected) , in which the contribu- 
tion of reflection is hardwired into the spectrum for a given 
geometry. However, in order to give a comparable number 
of degrees of freedom as in the CDID fits, we use xion to 
calculate the reflected emission only, and add a power law as 
the incident continuum so that the fraction Sl/2n is allowed 
to vary. We use the magnetic flares geometry to calculate 
the reflected emission from 10-100 Schwarzschild radii from 
a disc accreting at L/LEdd = 0.05 (see Fig. 1). The xion 
code includes relativistic smearing, but only for non-rotating 
black holes, so we switch this off and use the CONLINE convo- 
lution model as before. This gives a somewhat worse fit than 
before, with x 2 l v = 100/67. A better fit can be obtained by 
setting L/Z/Edd = 1, with x 2 l v = 86.1/67. This is due to the 
higher disc density expected at higher mass accretion rates. 
Both models give n n ~ 2.2 ±0.2 and q ~ 3.5 ±0.3 and again 
the best fit is reflection dominated. 

4.3 Inclination 

The inclination of the system is is not well known, but re- 
cent studies of the companion star indicate that its mini- 
mum value is 50 ± 3° (Orosz et al. 2004), showing that the 
assumption of 45° in the previous fits is probably too low. 

The filled squares in Fig. 3a show the increase in x 2 
with inclination for the fits for the CDID reflection models. 
Plainly the lower inclination fits are preferred by the data. 
However, we note that there is no minimum in x 2 a t 45° 
as seen from pexriv based reflection models (Miller et al. 
2002), so removing the only reason for using this inclina- 
tion. These fits become very unstable so we fix q — 3 and 
assume that the continuum is purely reflected emission (i.e. 
no intrinsic continuum visible) in order to produce the plots. 
However, we checked that these assumptions do not qualita- 
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Figure 2. MECS data fit to solar abundance reflection assuming an inclination of 45°. The (red in colour) dotted curve corresponds 
to the reflection, while the solid (cyan in colour) curve shows the sum of the reflection and the disc emission which is out of scale in 
these figures. Panel (a) shows the best fit using pexriv based reflection model, in which the reflector is highly ionized, and is extremely 
relativistically smeared. These relativistic effects strongly distort the intrinsic reflected spectrum. Panel (b) shows the very different 
reflected emission produced by CDID models, where the line and edge are intrinsically broad due to Comptonization. The best-fitting 
reflected spectrum is highly ionized, and extremely relativistically smeared. Panel (c) shows the same reflection model as in (b) but with 
no relativistic effects at all applied, so this fit is not dependent on inclination. The disc is again highly ionized, so the line is intrinsically 
broad and this intrinsic shape is a fairly good match to the data. This fit is only Ax 2 = 13 worse than the best fit (with the loss of 
2 degrees of freedom). Panel (d) shows the fit including ionized absorption from outflowing material. The inferred relativistic smearing 
reduces dramatically, so this fit is relatively independent of inclination 



tively affect the results, unless we allow for extremely cen- 
trally concentrated emission. We can obtain good fits with 
X 2 ~ 92 at higher inclinations up to 70°, but only with 
q ~ 15. 

At higher inclinations, the highly smeared solution be- 
comes a worse fit to the data because of the increasing 
blueshift of the relativistic transfer function due to the larger 
projected line of sight velocities (e.g. Fabian et al. 1989). 
The fit compensates for this by increasing the gravitational 
redshift i.e. by decreasing the disc inner radius and/or by in- 
creasing q. Since this figure is made for constant q = 3 then 
this requires that the inner radius decreases. However, this 
does not give enough gravitational redshift to compensate 
for the increased line broadening and most of the increase 
in x 2 comes from this increased blue extent of the spectral 
features. These extremely smeared fits then have a worse \ 2 
than the very different fits in which the line is mostly un- 
smeared, as described in the section above (see Fig. 2c). This 
solution is mostly independent of inclination as the intrinsic 
shape of the reflected emission below 10 keV is independent 
of inclination (Magdziarz & Zdziarski 1995), as are the rel- 
ativistic effects for large inner disc radii. Because these two 
solutions, the extremely smeared and almost unsmeared, are 
far apart in parameter space, they give a double minimum in 
X 2 space. We show the x 2 f° r the mostly unsmeared fits as 
the open triangles in Fig. 3(a), and it is plain that these are 
a much better description of the data for inclinations above 
45° than the highly smeared fits (filled squares). Fig. 3(b) 
shows the corresponding disc inner radius inferred for each 
solution as a function of inclination. Plainly, as the best fit 
switches between these two very different solutions, there is 
a corresponding dramatic jump in inferred disc inner radius. 



The residuals to the (mostly) unsmeared fit for an incli- 
nation of 45° in Fig. 2(c) show a clear feature at ~ 7.5 keV 
(the smeared fit in Fig. 2b also shows a similar feature). 
This looks like a fairly narrow absorption feature, and in- 
cluding a narrow, gaussian line (with best-fitting energy 
7.58 keV and equivalent width of ~100 eV) improves the 
fit by Ax 2 = 17.4, making it significantly better. The rel- 
ativistic smearing reduces dramatically, now being entirely 
consistent with q = 3, and with a larger disc inner radius, 
r- - 9+ 11 

' in — -J 4 • 



4.4 Narrow resonance line absorption 

Motivated by the fits above, we consider how such absorp- 
tion lines can occur. Observationally, strong resonance ab- 
sorption lines are often seen from high inclination Galactic 
binary systems, generally from iron Kq He- and H-like iron 
at 6.7 and 7.0 keV (e.g. Ueda et al 1998; Kotani et al. 2000; 
Lee et al. 2002; Ueda et al. 2004). However, if this absorbing 
material is outflowing then substantial blueshifts could be 
observed, and are seen in AGN (Pounds et al. 2003). If the 
feature at ~ 7.6 keV is an absorption line from iron then this 
implies a velocity of ~0.09c or 0.13c depending on whether 
the identification is with the H or He-like ion. 

We use the xstar photoionization models which are 
publicly available as xspec multiplicative model tables. We 
choose the ones calculated for power law illumination of a 
constant density photoionized column with turbulent veloc- 
ity of 100 km s' 1 . This gives an improved x 2 f° r incli- 
nations with the CDID reflection models (see Fig. 3c), and 
dramatically changes the inferred smearing of the reflected 
emission (Fig. 3d). The relativistic line is now always con- 
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Figure 3. Panel (a) shows the inclination dependance of x 2 f° r 
the MECS data fit with a disc blackbody and relativistically 
smeared CDID ionized reflection with q = 3. The best-fitting 
solution switches between a highly smeared model (filled squares, 
blue in colour) and a mostly unsmeared solution (open triangles, 
red in colour). Panel (b) shows the derived disc inner radius for 
each solution. For the extreme smearing, the increase in inclina- 
tion causes a large blueshift of the relativistic transfer function 
which is compensated by the fit by increasing the gravitational 
redshift, i.e. decreasing the inner radius. However, this also leads 
to increased broadening of the spectral features, and the total 
spectral shape does not match the data as well as at lower incli- 
nations. Conversely, for the mostly unsmeared fit, there is very 
little change with inclination in either x 2 or derived radius. Pan- 
els (c) and (d) show how these results change with the inclusion 
of outflowing ionized absorption (predominantly blueshifted El- 
and He-like iron resonance lines). This gives better fits at all in- 
clinations, for a disc inner radius which is always consistent with 
being larger than 6if, 9 . 



sistent with the illumination of q = 3 expected from gravi- 
tational energy release, and the inner disc radius is greater 
than QR g . 

The existence of several narrow features in the ionized 
absorber between 6.6-7 keV means that \ 2 space is complex, 
with multiple local minima. However, the column in all fits 
is ~ 10 22 cm" 2 of ionized material (£ ~ 100 erg cm s _1 ) 
outflowing at ~0.1-0.15c, with the best-fitting redshift of 
0.1511q'qq3. This model for an inclination of 45° is shown in 
Fig. 2(d). 

This effect is only seen when using the newer models 
of reflection from ionized material, which include the self- 
consistent Comptonization of the iron line and edge by the 
hot upper layers of the disc. Fig. 4 (grey, or green in colour, 
line) shows the best-fitting model using the pexriv reflected 
continuum and narrow Gaussian line, convolved with the 
LAOR relativistic profile. The sharp drop at ~7.2 keV seen 
in the data can be matched by the sharp blue wing of the 
relativistically smeared, intrinsically narrow line. By con- 
trast, the black line in Fig. 4 shows the best-fitting model 
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Figure 4. The grey (green in colour) line shows the best-fitting 
extreme smearing model where the intrinsic reflected emission is 
described by pexriv and a narrow gaussian line (Sec. 4.1). The 
sharp edge to the transfer function gives the steep drop of the 
blue wing of the iron line, while the deep edge in the ionised 
reflected continuum is smoothed into a broad feature. There is a 
fairly steep drop in the data at 7.2 keV so this model can match 
the curvature seen (see Fig. 2). The black line shows the best- 
fitting extremely smeared CDID model (Sec. 4.2). The line and 
edge are now mostly blended together, intrinsically broadened by 
Comptonisation in the intensely irradiated, hot upper layers of 
the reflecting slab. Relativistic broadening further smooths the 
spectrum, so there is no sharp drop at 7.5 keV. Since this drop 
is present in the data, this model is a worse fit despite it being 
physically more realistic. The arrow shows the energy inferred 
when fitting an additional absorption line to the data. 



using the extremely smeared, highly ionized CDID reflec- 
tion. The upper layer of the irradiated reflecting material 
is hot, so Comptonizes much of the line and edge, blending 
them together and making an intrinsically rather smooth 
spectrum (Ross et al. 1999). Convolution of this spectrum 
with relativistic effects makes it even broader, so there is 
no sharp change in curvature produced by the blue wing 
of the line. Since the data does have a sharp drop at 7.2- 
7.8 keV then these more physically self consistent models of 
reflection give a worse fit. The significant residuals in the 
spectrum at ~7.5 keV then are better fit by including out- 
flowing absorption, predominantly resonance iron K lines at 
rest frame energies of 6.7 and 6.95 keV. The drop seen in 
the data then no longer drives the relativistic line profile to 
extreme parameters and the derived smearing is completely 
compatible with a truncated disc. 



4.5 Summary of MECS fits 

There is a clear broad residual associated with iron in these 
BeppoSAX data from XTE J1650-500. Simple models of re- 
flection, which ignore Comptonization of the line, require 
that the relativistic effects are extreme, implying a spinning 
black hole and perhaps direct extraction of the spin energy. 
However, these models also require that the reflected emis- 
sion is highly ionized. This means that Comptonization of 
the line and edge is important, making these features in- 
trinsically broad and blended together. The strength and 
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significance of relativistic effects are then determined by the 
detailed modelling of the intrinsic reflected emission, which 
in turn depends on the (poorly known) vertical and radial 
structure of the disc. Nevertheless, the current best mod- 
els of ionized reflection always pick out parameters where 
the innermost radius of the emitting material is much closer 
than 6R g , and with radial emissivity rather steeper than ex- 
pected from gravitational energy release alone. However, the 
smearing derived from these (best available) reflection mod- 
els are not robust to changes in inclination, nor to the ad- 
dition of ionized absorption features. While either (or both) 
of these effects can resolve the discrepancy of the inner disc 
radius derived from the line shape with that inferred from 
the overall evolution of spectral and timing features, a fur- 
ther conflict remains which is that the fits favour a reflection 
dominated solution. 



5 MECS AND PDS 2-200 KEV: COMPLEX 
CONTINUUM MODELS 

5.1 Simple reflection models 

To get more constraints on the amount of reflection requires 
data at higher energies, so we include the simultaneous PDS 
data to extend the spectrum up to 200 keV. However, the 
spectrum is now much more sensitive to details of the contin- 
uum modelling due to the wide bandpass, and it is obvious 
that a power law can no longer be used. Even an expo- 
nentially cutoff power law is not a good approximation to 
the shape of a thermally Comptonized spectrum over such 
a wide bandpass. The rollover at energies close to the elec- 
tron temperature, kT e , is rather sharper than an exponen- 
tial (Gierlihski et al. 1997). Also, there are indications that 
the spectra can be more complex than predicted by just a 
single temperature thermal Compton component, with sig- 
natures of both thermal and non-thermal electrons even in 
the low/hard state (McConnell et al. 2000; Ibragimov et al. 
2005). 

We use the eqpair continuum model (Coppi 1999) in 
order to model the curvature of the Comptonized contin- 
uum. This computes the Comptonization of seed photons 
by a steady state electron distribution which is calculated 
self-consistently from balancing an initial heating rate with 
Compton and Coulomb cooling processes. The code allows 
the input heating to be either non-thermal, with electrons 
injected with a power-law distribution of Q(-~i) oc 7 _s , be- 
tween Lorentz factors 7 m in and 7 max , or thermal, or some 
combination of both. However, even if all the electron heat- 
ing is non-thermal, the self-consistent electron distribution 
(and hence output Comptonized spectrum) can be mainly 
thermal with only a small non-thermal tail if Coulomb colli- 
sions dominate the electron cooling. This happens where the 
power injected in the hard electrons, lu (parameterized as a 
compactness I = Lot / Rm e c 3 ; see e.g. Coppi 1999) is larger 
than that injected in soft photons, £ s . Conversely, where the 
soft photons dominate then Compton cooling is most impor- 
tant, and the electron distribution retains a mainly power- 
law shape with only a small thermal component (Poutanen 

6 Coppi 1998). 

As well as calculating the curvature of the Comptonized 
spectrum at high energies, eqpair also computes the low- 
energy cutoff due to the seed photon distribution. This cutoff 



is typically at a few kT BCO d , so is within the observed band- 
pass for the MECS data. This curvature can make a large 
difference to the inferred disc temperature and luminosity 
(see e.g. Done, Zycki & Smith 2002), and could potentially 
distort the inferred spectrum in the red wing of the iron 
line. We assume that the seed photons are the same as the 
observed disc photons in all the fits below, and fix £ s — 10. 

This continuum complexity means that we cannot eas- 
ily use the CDID reflection models since these are tabulated 
assuming a cutoff in the continuum at a fixed energy of 
100 keV. The assumed high energy continuum shape affects 
the reflected emission at energies above ~ 15 keV due to the 
predominance of Compton down-scattering. Instead we use 
the reflection model in eqpair (which is based on PEXRIV 
but reflecting the self-consistent continuum from eqpair), 
together with a gaussian line. We relativistically smear both 
reflected continuum and line with the LAOR kernel, now ex- 
tending the bandpass for the convolution calculation from 
0.2-1000 keV, and assume an inclination of 45°. Purely ther- 
mal electron heating gives an (almost) acceptable fit, with 
xl = 99.9/79 for a reflection dominated (ft/27r = 2.4±J;|), 
highly smeared fit (q = 3.8 ± 0.2, r in = l-8±g;f ). How- 
ever, purely non-thermal injection models give a significantly 
better fit, with \ 2 = 82.7/78 (there is one extra free pa- 
rameter which is the spectral index of the power law elec- 
tron injection). The subtle difference in continuum curvature 
means that the solution is no longer reflection dominated, 
with Q./2-k = 0.72l°'Jg though it is still strongly smeared 
(q = 3.5± | and r in = 2.2±g;|). The MECS and PDS spec- 
trum and residuals for this best-fitting non-thermal injection 
model are shown in Fig. 5(a). 

5.2 Complex reflection models 

Section 4.2 showed that there are large differences in the 
model spectra for ionized reflection. While the CDID reflec- 
tion models concentrate on the spectra below ~12 keV, com- 
puting the effects of photoionization and Compton up and 
down-scattering at low energies, they do not accurately treat 
Compton scattering at higher energies, nor include angle de- 
pendent effects, nor allow for different illuminating continua. 
This contrasts with the pexriv models of reflection which 
have very approximate ionization balance (Done et al. 1992) 
and neglect Comptonization below 12 keV, but treat Comp- 
tonization at higher energies very accurately, and give the 
reflected emission as a function of observed inclination an- 
gle, pexriv is based on the convolution kernels given by 
Magdziarz & Zdziarski (1995), so can be extended to work 
with any illuminating continuum. 

Plainly what is needed is a combination of these models, 
using the better CDID photoionization/Compton scattering 
of reflection at low energies and the pexriv accurate Comp- 
ton down-scattering at higher energies. The most flexible 
way to incorporate this is as a convolution model in xspec, 
so that the reflected spectrum can be calculated for any arbi- 
trary continuum. We have coded this by finding the pexriv 
cross-section at 10 keV which gives the same 12-14 keV re- 
flected spectrum gradient as the CDID models for a given 
ionization parameter and illuminating spectral index. The 
corresponding pexriv Greens functions are used to form the 
spectrum above 12 keV. Below 12 keV the Greens functions 
are assumed to collapse to delta functions, given by the ratio 
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of reflected to incident flux in the CDID models, rescaled to 
match the normalization of pexriv at 14 keV. This model 
is hereafter called refbal. 

We use this together with the eqpair continuum model, 
and get a somewhat worse fit than for the pexriv models, 
with x 2 = 102.5/80 for an illuminating continuum from non- 
thermal injection. This is shown in Fig. 5(b) and by contrast 
with the pexriv reflection, it requires a larger fractional 
contribution from reflection, with 0,/2-k = 2.01J; 2 . However, 
both pexriv and CDID based models require that reflection 
is highly smeared, with n n = 2-3. 

We then include the narrow absorption features from 
XSTAR ionized absorption. As with the MECS data alone, 
this dramatically reduces the derived relativistic smearing 
to being consistent with a disc truncated above r ln = 6 
for the CDID reflection models. However, the fit still has 
0./2n = 2.41q 4 i.e. is reflection dominated, although there 
is a significant contribution from the intrinsic continuum 
component as well. 

5.3 Summary of MECS and PDS fits 

The inclusion of the PDS data up to 200 keV plainly gives 
more constraints on the amount of reflection than the MECS 
data alone. With data only up to 10 keV the amount of re- 
flection is very poorly defined, and the spectrum can be pure 
reflected emission with no intrinsic continuum. Including the 
higher energy data shows that there is most probably true 
continuum present, though the derived reflected fraction is 
sensitive to details of how the reflected and continuum emis- 
sion are modelled. With the pexriv based reflection mod- 
els, the data are consistent with fl/2n < 1 as expected for a 
truncated disc, while the more physical CDID based reflec- 
tion models have substantially more reflection present, with 
Q./2-k ~ 2. Both these require extreme smearing and emis- 
sivity, but like the MECS-only fits, the smearing is again 
sensitive to the inclusion of narrow absorption lines from an 
ionized outflow for the CDID reflection models. 



6 DISCUSSION 

The bright low/hard state spectrum seen by BeppoSAX 
from the transient black hole XTE J1650-500 is certainly 
consistent with extremely smeared, reflection dominated 
emission as claimed by Miniutti et al. (2004). However, it 
is also consistent with much less smearing if there is some 
ionized, outflowing absorption which gives narrow resonance 
lines at ~7.5 keV. Reducing the amount of smearing makes 
the reflected iron line consistent with the expectations of 
truncated disc models, removing the conflict between the 
position of the inner edge of the cool disc derived from 
the broad-band spectral shape and QPOs. This is directly 
testable with high resolution data around the iron line. 

Such absorbing material is often seen in X-ray binaries, 
especially those at fairly high inclination angles (Ueda et 
al. 1998; Kotani et al. 2000; Lee et al. 2002; Ueda et al. 
2004), including XTE J1650-500 (Miller et al. 2004a). The 
one difference here is that the outflowing material has veloc- 
ity ~4.5xl0 4 km s _1 , much higher than the ~500 km s _1 
typically seen in these objects. A low velocity wind would 
be launched from fairly large radii in the disc, as expected 



for line driving from the UV emitting region of the GBH 
disc (Proga & Kallman 2002). By contrast, the much higher 
velocities required here would need to be associated with a 
launch from the inner accretion disc, where line driving is 
not effective due to the high ionization. This suggests rather 
that the wind would have to be magnetically driven, and/or 
be associated with the jet. We note that these data are taken 
shortly after the time of the maximum in radio emission 
(Corbel et al. 2004), so the inferred very fast outflow could 
be associated with the jet ejection. 

We can estimate the mass outflow rate. For a constant- 
velocity outflow the mass outflow rate at a radius 7? is 

Mout = iirR 2 nm p fv, (1) 

where n is the number density of the outflowing material, 
v is its velocity, / = Q./An is the fractional solid angle in 
which the outflow is directed and m p is the proton mass. 
In a steady state M ou t is constant for all radii. The column 
density measured by the observer is 

N H = [°°ndR= A 1 ^° ut , (2) 
J ATYfrrip R in v 

therefore 

Mout = 4-wfm p N H R in v. (3) 

For N H = 10 22 cm" 2 , R in = WR g , v = 0.15c and / < 
1 we find M out < 1.4 x 10 16 g s _1 « 10~ 3 Af Edd (assuming 
accretion efficiency of 0.1 and black hole mass of 10 Mq). 
This is only about 1 per cent or less of estimated inflow mass 
accretion rate and constitutes a very feasible number. 

Our spectral fitting also shows that this absorption does 
not make the data consistent with all the expectations of a 
truncated disc. The best current models of reflection based 
on irradiation of a slab of constant density give reflection 
fractions which are much larger than expected for a trun- 
cated disc geometry. This is not necessarily the case using 
the older pexriv reflection models, which can give much 
smaller reflection fractions. The difference is due to Comp- 
tonization of the line and edge in the reflected emission as 
it travels through the upper layers of the disc (Ross et al. 
1999). This effect is neglected in pexriv, so the iron edge is 
rather sharp. On the other hand, the CDID reflection is not 
calculated for the disc temperatures appropriate for Galac- 
tic binaries i.e. 0.5-1 keV. In such hot material, collisional 
processes alone will ionize iron up to the K shell (Davis et 
al. 2005). By contrast, the CDID models do not include col- 
lisional ionization, and thus require intense illumination to 
ionize iron up to the K shell to match the data. This pulls 
the temperature of the upper layers to the Compton temper- 
ature of 1.5-2 keV, for spectra with F ~ 1.9-2 (Ross et al. 
1999; Nayakshin et al. 2000). This is a factor 2-3 larger than 
the observed disc temperature in these data. Thus while the 
pexriv models underestimate the width of the line and edge 
by having no Compton scattering below 12 keV, the CDID 
models may overestimate it for Galactic binaries. From our 
spectral fits, the derived reflected fraction is sensitive to the 
detailed shape of this edge. Better models of reflection for 
conditions appropriate to Galactic binary discs will enable 
this effect to be quantified, removing this as a source of un- 
certainty. 

The amount of reflection is also sensitive to the detailed 
form assumed for the continuum. This is not easily defined 
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Figure 5. MECS and PDS data fit with a disc blackbody (dashed curve, green in colour) and EQPAIR continuum (long-dashed, blue in 
colour), together with the reflection (dotted, red in colour) from the disc. Panel (a) shows the model with simple reflection, panel (b) 
with CDID reflection and panel (c) the same, but with ionized blueshifted absorption added. 



even with broad bandpass data. Obviously, the wider the 
bandpass, the more constraints there are on the form of 
the continuum curvature, but also the more the spectrum is 
sensitive to the detailed form of continuum assumed. Broad 
bandpass data suggest that the spectrum has convex cur- 
vature which can be modelled by an additional soft excess 
component (e.g. Di Salvo et al. 2001; Ibragimov et al. 2005). 
This makes the intrinsic continuum at higher energies rather 
harder than expected from the observed 2-10 keV emis- 
sion, reducing the amount of reflection required to match 
the data at high energies. Including such a component in 
the model makes the reflected fraction almost completely 
undetermined and consistent with a truncated disc with 
the refbal and ionized absorber models. Such additional 
continuum complexity is not required statistically from the 
spectral data alone (even with the moderate spectral res- 
olution and broad bandpass of BeppoSAX). However, en- 
ergy resolved variability could give additional constraints as 
such data can fairly unambiguously pinpoint the seed pho- 
ton energy seen by the Comptonizing plasma (Gierlihski & 
Zdziarski 2005), so could disentangle the shape of the soft 
emission. 



7 OTHER EXTREME LINE OBJECTS 

Most detections of extreme broad iron lines in GBH are in 
the very high state (alternatively named the steep power law 
state) rather than in the low/hard state examined here (see 
e.g. Miller et al. 2004a). This (typically rather high mass 
accretion rate) spectral state shows strong disc emission as 
well as a strong, steep high-energy tail, while the power spec- 
trum of the variability shows a strong QPO which is close to 
its maximum frequency (e.g. McClintock & Remillard 2006). 
Thus both spectral and timing properties imply that the in- 
ner disc does extend close to the black hole (quite how close 
is difficult to determine due to the complexity of the spectra: 
Kubota & Done 2004; Done & Kubota 2005). Thus there are 



no obvious conflicts as for the low/hard state. There are ob- 
jects where the disc dominated spectra imply moderate spin 
(from the Ldisc oc T^ax relation) while the iron line smear- 
ing implies extreme spin but these could be explained by the 
accretion flow extending below J? ms in the very high state 
(Reynolds & Begelman 1997). 

However, there is a more subtle conflict with the ex- 
treme line in these data which is that the very high state 
spectra typically show strongly Comptonized emission. The 
broad band spectral shape seen in e.g. the simultaneous 
ASCA-RXTE-OSSE data from XTE J1550-564 implies that 
the inner disc is almost completely covered by material with 
optical depth 2-3 (Done & Kubota 2005; Gierlihski & Done 
2003). Thus the spectral shape implies that very little of the 
inner disc is seen directly, so it cannot make much contribu- 
tion to the line emission. Yet the same ASCA data require an 
extreme broad line when fit by pexriv based models (Miller 
et al. 2004b). Again, we suggest that resonance line absorp- 
tion in an outflowing wind can modify the derived line width 
and remove these conflicts in the very high state spectra. To 
test this, we fit the publicly available standard products Gas 
Imaging Spectrometer data, as this does not suffer from the 
pileup problems of the Solid-State Imaging Spectrometer for 
such a bright source. The relativistic smearing required for 
a pexriv reflected continuum and narrow line are extreme, 
with r in = 1.6±° ;l and q = 3.7 ± 0.3, with X 2 = 796.5/742, 
similar to the parameters found by Miller et al. (2004b). 
Replacing this with CDID reflection gives similar extreme 
relativistic smearing but a worse fit (x 2 = 824.6/744). In- 
cluding XSTAR ionized absorption improves the fit {\ 2 = 
805.4/741), and as in XTE J1650-500 this reduces the rel- 
ativistic effects to the point where ri„ = 31_is (no upper 
limit) with q fixed at 3. 

The extreme broad lines inferred in some AGN are sim- 
ilarly potentially affected by absorption complexity chang- 
ing the line parameters. Here there are as yet no obvious 
problems with this interpretation, especially as extreme spin 
might be expected given the hierarchical growth of black 
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holes through accretion (Volonteri et al. 2005). Nonethe- 
less, if the apparently extreme GBH iron lines are actually 
consistent with low spin, normal illumination rather than 
extraction of rotational energy from an extreme spin black 
hole, then it is likely that the AGN discs behave similarly. 
We note that the recent Chandra spectrum of the most se- 
cure extreme relativistic line, that in MCG-6-30-15, reveals 
the presence of ionized absorption at a similar column and 
ionization as that required here to reduce the relativistic 
smearing to levels consistent with a truncated disc (Young 
et al. 2005). Young et al. (2005) checked that this absorption 
did not change the extreme smearing deduced from a LAOR 
description of the line, but this is unsurprising due to the 
unphysically sharp blue wing of the line produced in such 
models (see Fig. 4). Our fitting showed that the relativis- 
tic line parameters only changed significantly when using 
the CDID highly ionized reflection models as these have a 
much smoother intrinsic spectrum around the iron line/edge 
due to Comptonization. We urge testing the robustness of 
the extreme relativistic effects to the inclusion of ionized 
absorption with more physical ionized reflection models. 



ate columns of highly ionized absorbing material (predom- 
inantly resonance lines of H- and He-like iron) can always 
reduce the inferred relativistic smearing of a highly ionized 
reflected spectrum. This could then remove the more sub- 
tle conflict seen in the very high state GBH spectra, where 
the disc should be down close to the last stable orbit (so 
strong relativistic smearing is expected) but the spectral 
shape shows that the inner disc is mostly covered by an 
optically thick, hot corona. We stress that the presence of 
absorption at this level is directly testable with high resolu- 
tion data. 

While there are no such direct conflicts for AGN, we 
note that the most secure extreme line candidate, MCG-6- 
30-15, has highly ionized reflection (e.g. Wilms et al. 2002) 
and a recent Chandra HETG spectrum shows absorption 
at the same column and ionization (Young et al. 2005) as 
inferred here for XTE J1650-500. We urge a re-analysis 
of these data to see whether the extreme line parameters 
are robust to the presence of absorption when using the 
more physical reflection models which include self-consistent 
Comptonization of the line and edge as opposed to models 
in which the line is intrinsically narrow. 



8 CONCLUSIONS 

The truncated disc/hot inner flow geometry is very success- 
ful in qualitatively (and sometimes quantitatively) explain- 
ing many disparate properties of the low/hard state in GHB. 
If the radius at which the hot flow evaporates from cool disc 
progressively decreases as the source approaches the transi- 
tion to the soft state then this can account for the observed 
evolution in spectral and timing properties. However, if there 
really are detections of substantial amounts of iron line emis- 
sion from cool disc-like material down at the last stable orbit 
in the low/hard state, then the truncated disc idea is simply 
wrong. Some of the extreme line emission might be explained 
by clumps from a disrupted disc embedded in the hot flow, 
but this would give rise to only a small amount of line since 
the covering fraction of the clumps would be rather less than 
from a disc, and these are embedded in a hot flow with opti- 
cal depth of around unity, so Compton scattering suppresses 
the line. 

We show that the bright low/hard state spectrum seen 
by BeppoSAX from the transient black hole XTE J1650- 
500 can indeed be interpreted as requiring strong, extremely 
smeared reflected emission from a disc-like structure extend- 
ing down to ~ 2R g as shown by Miniutti et al. (2004). 
These data clearly illustrate the conflict between this and 
the position of the inner edge of the cool disc derived from 
the broad band spectral shape and QPOs. The small inner 
radius inferred from the line width also conflicts (though 
more weakly) with the moderate spin inferred both from 
theoretical models of stellar collapse and from the observed 
Ldisc oc T^ ax relation in this object. We show that all these 
conflicts can be removed if relativistic effects on the line and 
reflected continuum are overestimated due to the presence 
of absorption line features in the spectrum. The derived line 
width is then easily consistent with normal emissivity and 
a disc truncated at ~10-20 R g (so giving no information on 
the black hole spin). The reflection feature are most prob- 
ably shaped by strong, but not necessarily extreme, grav- 
ity. We suggest that this effect is robust, i.e. that moder- 
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